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From  mid-January  1974  through  March,  two  sets  of  LASA  A-D  ring  beams  were 
ormed  in  parallel  by  the  on-line  LASA  Detection  Processor  (DP  at  SDAC) . One 
set  contained  299  LASA  array  beams  using  phased  subarray  sums;  the  other  con- 
tained array  beams  to  the  same  299  regions  using  unphased  subarray  sums.  In 
six  experiments  we  used  seismograms  recorded  by  the  inner  1,  4,  7,  10,  13 
and  all  16  sensors  per  subarray  to  form  the  unphased  sums  for  the  LASA  beams. 
Each  of  the  LASA  beams  so  constituted  was  compared  to  the  LASA  beam  containing 
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phased  sums  of  16  traces  per  subarray.  The  objective  of  these  experiments  was 
to  determine  the  effect  on  array  detection  performance  of  using  unphased  sums. 

A comparison  of  the  signal-to-noise  (S/N)  ratios  of  the  on-line  detections 
from  the  parallel  beam  sets  shows  that  the  average  S/N  loss  on  the  beams  con- 
taining unphased  subarray  sum  traces  varies  from  2.4  db  for  16  channels  to 
i.6  db  for  one  channel.  This  suggests  a change  in  threshold  varying  between 
.1  and  0.3  magnitude  units;  and  the  relative  numbers  o£  events  detected  by 
the  phased  and  unphased  subarray  beams  support  this  interpretation.  This  is, 
however,  an  average  over  the  teleseismic  distance  range;  the  loss  is  greater 
for  ciose-in  distances  than  for  large  distances  and  is  partly  due  to  higher 
noise  levels  because  of  the  smaller  number  of  sensors  in  the  smaller  subarray, 
separation  and  discussion  of  these  two  effects  as  a function  of  frequency 
leads  to  the  conclusion  that  10  sensors  per  subarray  strikes  a satisfactory 
balance  between  high-frequency  signal  loss  and  reduction  of  the  detection 
threshold.  There  is,  however,  no  sharp  change  at  any  particular  number  of 
elements,  and  other  choices  are  possible. 
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ABSTRACT 


From  mid- January  1974  through  March,  two  sets  of  LASA  A-D  ring  beams  were 
formed  in  parallel  by  the  on-line  LASA  Detection  Processor  (DP  at  SDAC) . One 
set  contained  299  LASA  array  beams  using  phased  subarray  sums;  the  other  con- 
tained array  beams  to  the  same  299  regions  using  unphased  subarray  sums.  In 
six  experiments  we  used  seismograms  recorded  by  the  inner  1,  4,  7,  30,  13 
and  all  16  sensors  per  subarray  to  form  the  unphased  sums  for  the  LASA  beams. 
Each  of  the  LASA  beams  so  constituted  was  compared  to  the  LASA  beam  containing 
phased  sums  of  16  traces  per  subarray.  The  objective  of  these  experiments  was 
to  determine  the  effect  on  array  detection  performance  of  using  unphased  sums. 

A comparison  of  the  signal-to-noise  (S/N)  ratios  of  the  on-line  detections 
from  the  parallel  beam  sets  shows  that  the  average  S/N  loss  on  the  beams  con- 
taining unphased  subarray  sum  traces  varies  from  2.4  db  for  16  channels  to 
5.6  db  for  one  channel.  This  suggests  a change  in  threshold  varying  between 
0.1  and  0.3  magnitude  units;  and  the  relative  numbers  of  events  detected  by 
the  phased  and  unphased  subarrcy  beams  support  this  interpretation.  This  is, 
however,  an  average  over  the  teleseismic  distance  range;  the  loss  is  greater 
for  close-in  distances  than  for  large  distances  and  is  partly  due  to  higher 
noise  levels  because  of  the  smaller  number  of  sensors  in  the  smaller  subarrav. 
Separation  and  discussion  of  these  two  effects  as  a function  of  frequency 
leads  to  the  conclusion  that  10  sensors  per  subarray  strikes  a satisfactory 
balance  between  high-frequency  signal  loss  and  reduction  of  the  detection 
threshold.  There  is,  however,  no  sharp  change  at  any  particular  number  of 
elements,  and  other  choices  are  possible. 
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INTRODUCTION 


Tht  objective  of  this  study  is  to  determine  the  effects  of  using  unphased 
subarrav  sums  on  the  detection  and  discrimination  performance  of  the  Large 
Aperture  Seismic  Array  (LASA)  in  Montana  and  to  make  preliminary  recommenda- 
tions of  the  number  of  sensors  to  use  to  form  the  unphased  sums. 

It  is  anticipated  that  the  50  kilobit  line  capability  from  Montana  to 
the  SDAC  will  be  reduced  to  a 4.8  kilobit  line  in  the  near  future.  To  conform 
to  this  reduced  transmission  rate,  the  data  output  from  LASA  can  be  decreased 
by  transmitting  unphased  subarray  sums  from  Montana  from  each  of  the  thirteen 
remaining  subarrays  instead  of  transmitting  the  individual  outputs  from  the 
208  sensors  still  in  operation. 

As  programmed  for  this  experiment,  the  detection  process  in  Partition  I 
serves  as  a based  of  comparison  for  the  unphased  sums  (calculated  in 
Partition  II).  In  Partition  I we  form  five  16-sensor  subarray  phased  sums 
equi-spaced  in  azimuth  at  a velocity  of  15.5  km/sec  (0.0645  sec/km)  for  each 
of  the  thirteen  subarrays.  The  subarray  beams  are  combined  into  the  full 
array  beamset  LBS  140  which  selects  from  each  subarray  the  beam  which  is 
closest  in  velocity  space  to  the  beaming  point  of  the  full  array  beam.  The 
299  beaming  points  in  the  full  array  beam  are  hexagonally  packed  in  inverse 
velocity  space  with  a spacing  of  0.001316  sec/km  out  to  a radius  of  0.132 
sec/km  or  7.6  km/sec. 

These  beams  are  then  examined  by  a detection  process  described  by 
Chang  (1974).  For  ease  of  reference,  we  shall  brief 1>  review  the  process. 

The  subarrav  beams  in  this  experiment  have  been  filtered  0.8-2. 5 Hz 
instead  of  0.9-1. 4 Hz  as  was  previously  implemented  in  the  SDAC  system. 

Bacause  of  the  broader  pass-band,  pulse-like  signal  wave-forms  are  not  so 
extended  in  time  after  passing  through  the  filter. 

This  fact  suggests  that  changes  in  other  detection  parameters  might  be 
desirable.  Such  changes  were  made  and  are  discussed  below.  The  first  step 
in  detection  is  to  accumulate  a short-term  average  (STA)  of  the  noise  by 
averaging  over  1.8  seconds  of  the  rectified  beam  of  interest.  Successive 
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overlapping  STA  values  are  computed  every  0.6  seconds.  From  these  STA's  a 
long-term  average  is  computed  every  third  step,  or  once  in  1.8  seconds,  by 
exponentially  weighting  the  previous  LTA  value  and  adding  the  current  STA 
value.  The  detection  algorithm  tests  if  the  STA/LTA  ratio  is  greater  than 
a S/N  acceptance  threshold  (10  db  in  these  experiments)  each  time  an  STA  is 
computed. 

The"  are  two  parts  to  the  detection  algorithm  for  LASA.  The  first  part 
is  the  signal  to  noise  ratio  (S/N)  threshold  teSL  which  determines  the  size 
and  duration  of  the  signal.  When  the  ratio  of  the  STA/lta  exceeds  the  fixed 
threshold  value  of  10  db  for  the  duration  of  Q out  of  Q'  (2/2  in  the  experi- 
ment; changed  from  3/3  as  previously  implemented)  consecutive  tests,  the 
signal  arrival  is  declared  "on"  on  that  beam.  After  the  beam  is  turned  on, 
the  end  of  the  signal  arrival  is  declared  when  the  S/N  ratio  of  the  beam 
becomes  lower  than  the  turn-off  threshold  of  7 db.  The  LTA  computation  is 
stopped  when  the  beam  is  "on".  Note,  however,  that  the  LTA  may  be  contami- 
nated by  the  signal  because  if  the  LTA  is  being  updated  at  the  first  success- 
ful threshold  test  (one  chance  in  three),  then  the  STA  being  added  to  the  LTA 
will  contain  some  signal.  (While  this  should  not  affect  detection  thresholds 
since  weak  events  cannot  significantly  affect  the  LTA,  it  will  lower  the 
reported  S/N  values  for  large  events  one-third  of  the  time.) 

The  second  part  of  the  detection  algorithm  is  the  spatial  consistency 
test.  This  determines  the  consistency  of  the  seismic  signal  in  both  azimuth 
and  velocity  by  seeking  the  maximum  STA  beam  and  checking  if  the  previous 
maximum  was  found  within  the  distance  of  AU  beams  (AU  = 2 beams)  from  the 
current  maximum  beam.  When  this  condition  is  satisfied  for  ? (P  = 3 in  this 
experiment;  changed  from  4 as  previously  implemented)  consecutive  times,  the 
signal  arrival  is  declared  on  the  beam  with  the  highest  STA  value  during 
these  P consecutive  tests.  If  that  beam  is  "on"  by  virtue  of  having  passed 
the  threshold  test  then  both  parts  of  the  detection  algorithm  have  been 
passed  and  a detection  is  declared. 

To  carry  out  our  experiments,  this  detection  system  was  placed  in 
Partition  I of  the  SDAC  system.  In  Partition  II  we  placed  an  identical  system 
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except  that  only  one  beam,  an  infinite  velocity  beam,  was  formed  at  the 
subarray  level;  and  a variable  number  of  sensors  were  summed  into  the  subarray 

beam  depending  on  the  experiment  under  consideracion  as  illustrated  in 
Table  I. 

The  complete  printed  Detection  Processor  (DP)  outputs  are  available  for 
both  partitions;  and  we  have  compared  the  number  of  events  which  crossed  the 
12  db  threshold  and  were  confirmed  by  the  analysts  in  the  two  partitions. 

As  it  happened  there  was  never  a case  of  a detection  above  12  db  in  Parti- 
tion II  and  not  in  Partition  I which  was  accepted  as  an  event  by  the  analysts. 
Thus  we  have  tabulated  the  percentage  of  analyst-accepted  Partition  I events 
above  12  db  which  were  also  detected  in  Partition  II  above  12  db.  We  have 
also  tabulated  the  difference  in  (S/N)  for  all  accepted  events  with  (S/N) 
greater  than  10  db  in  Partition  II.  This  is  as  low  as  it  is  possible  to 
tabulate  because  lower  db  values  are  not  printed  out. 

No  comparison  was  made  of  the  false  alarm  rate  between  the  two  partitions; 
however,  there  is  every  reason  to  expect  that  they  would  be  identical.  By 
reference  to  Blandford  and  Clark  (1971)  the  intersensor  noise  correlation 
inside  a subarray  is  very  low  in  the  0. 8-2.0  Hz  pass-band  for  both  infinite 
velocity  and  12  km/sec  beams.  The  noise  correlation  between  subarrays  is 
generally  regarded  as  zero,  and  in  general  one  does  find  */N  noise  reduction 
by. beamforming  subarrays,  Hartenberger  and  van  Nostrand  (1970).  Thus,  no 
matter  what  the  subarray  beam,  and  no  matter  how  many  sensors  in  the  subarray, 
one  expects  the  noise  processes  to  be  identical  among  partitions  and  among 
experiments  except  for  the  mean-square  level.  The  mean-square  level  will, 
however,  not  affect  tne  false  alarm  rate  since  the  threshold  is  a relative 
one;  comparing  a short-term  average  of  the  noise  to  a long-term  average. 
Examination  of  a day’s  DP  output:  February  20,  1974  showed  no  difference 

between  partitions  in  the  number  of  threshold  crossings  flagged  by  the  event 
processor  as  sidelobe  detections  or  as  later  phase  arrivals.  In  any  event, 
such  an  increase  would  result  in  no  increased  load  on  the  analysts. 

It  would  have  been  possible  to  run  each  experiment  on  the  same  data  by 
use  of  the  off-line  DP  system.  If  this  had  been  done,  a more  controlled 
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comparison  would  have  been  possible  between  experiments.  Unfortunately  the 
otf-line  DP  runs  at  one-half  real  time;  and  the  360/40B  required  for  the  off- 
line runs  is  available  only  6 hours  per  day  on  the  average.  Thus  to  run  the 
same  number  of  experiments  would  have  required  eight  times  as  long;  or  2-3 
years  instead  of  3-4  months,  see  Table  I.  Off-line  DP  tests  also  encounter 
severe  tape  reading  problems,  Chiburis  et  al.  (1974). 

The  required  results  could  also  have  been  obtained  at  very  low  conq u ra- 
tion cost  simply  by  beamforming  a few  large  events  for  each  distance  interval 
and  region  of  interest.  Experience  has  shown  that  the  standard  deviation  of 
signal- to-noise  improvement  estimates  obtained  by  this  method  is  very  low; 
and  that  only  a few  events  are  needed  to  define  the  mean.  (In  contrast,  of 
course,  a very  large  number  of  events  are  required  to  accurately  determine 
a threshold  by  the  technique  of  plotting  probability  of  detection  or  number 
of  detections  versus  magnitude.) 

However,  it  is  desired  to  test  the  actual  system  in  operation  so  that 

any  unsuspected  effects  of  the  complicated  detection  algorithms  just  discussed 
would  be  revealed. 


-4- 


THEORETICAL  CHANGES 


In  Table  II  we  see  some  of  the  principal  results  of  experiments  1-6. 

We  see  in  Table  I that  the  subarrays  were  filled  from  the  center  out  as  the 
number  of  the  experiment  increased.  It  would  seem  that  these  are  the  sub- 
array designs  of  interest  since  in  this  way  the  signal  loss  on  unphased  sums 
due  to  signal  misalignment,  and  the  costs  of  cable  and  instrument  maintenance 
increase  monotonically . Also,  as  we  shall  see,  the  most  closely  spaced 
sensors  contribute  in  all  subarray  designs  to  reduction  of  noise.  (It  is  of 
interest  in  this  regard  that  the  original  array  design  of  the  Geneva  experts 
called  for  an  array  of  3 km  diameter  because  zero-delay  analog  summation  of 
the  individual  sensors  was  envisaged  as  the  processing  scheme  of  choice.) 

In  each  of  the  experiments  the  average  distance  to  the  events  lay 
between  63  -70°  with  a mean  of  67°;  and  the  average  waveform  period  in  each 
experiment  was  0.9  seconds.  To  the  extent  that  the  stability  of  these  sta- 
tistics indicates  that  the  sample  populations  range  of  113-379  events  can 
be  taken  as  a satisfactory  sample  of  the  population  of  interest,  we  may 
obtain  a crude  estimate  of  the  threshold  change  from  the  percent  detection 
in  the  second  partition. 

It  is  generally  assumed  that  the  probability  of  detection  is  a function 
cf  the  signal-to-noise  ratio,  that  is  to  say  a function  of  the  difference  in 
the  logarithm  of  the  signal  and  noise.  Blandford  and  Wirth  (1973)  have 
shown  that  this  function  may  be  approximated  for  automatic  detectors  by  a 
step  from  0 to  1.0  proability  at  specified  (S/N).  Then,  for  any  particular 
array,  the  probability  of  detection  is  equal  to  the  probability  of  attaining 
a (S/N)  value  greater  than  a particular  threshold  value.  For  a fixed  LASA 
magnitude  at  a fixed  distance,  the  magnitude  is  proportional  to  the  log 
amplitude.  Since  the  mean  noise  amplitudes  have  been  found  by  many  workers 
to  be  distributed  log-normally,  the  probability  of  detection  of  a fixed 
is  given  by 

P(mb”lJ)  = P(fmb"b]/2o)  = <f(  [m^-p  ]/2o) 

where  u = B(A)  + log  a + log  r combines  the  logarithm  of  the  (S/N)  threshold 
value  r;  the  logarithm  of  a,  the  median  noise  level;  and  B(A)  the  distance 
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amplitude  factor  for  P waves;  o is  the  standard  deviation  of  the  log  noise 
and  <p  is  the  cumulative  normal  distribution  function.  If  we  neglect  the 
variation  of  log  amplitude  with  distance  for  fixed  for  the  events  in  the 

LASA  detection  list  we  may  use  the  above  formula  for  all  events. 

Then  if  the  number  of  events  with  magnitude  is  given  by 

-brn 

n = C 10  b 

then  the  total  number  of  detected  events  is 

a -bim 

N = C / p(m^-n)  10  bdrm 

— 00 

and  a transformation  of  the  variable  of  integration  q = nt  - p shows  that: 

-b(p  -p  ) 

n(p1)/n(p2)  =10  1 2 . 


Thus  if  we  know  b,  then  from  the  change  in  percent  events  detected  in 
Partition  II  compared  to  those  in  Partition  I,  we  can  determine  the  change 
in  the  "effective"  threshold. 

From  Figure  8 in  Chang  (1974),  the  cumulative  recurrence  curve  for 
1500  LASA  EP  amplitudes  we  obtain  b = 1.06;  from  Dean  (1971)  we  obtain 
b = 0.9  from  a recurrence  curve  of  500  LASA  values.  Let  us  choose 
b = 1.0.  By  application  of  the  above  analysis,  we  arrive  at  column  d in 
Table  II. 

To  determine  the  (S/N)  loss  in  each  experiment  we  averaged  the  difference 
in  log  (S/N)  in  each  partition  for  all  events  detected  in  Partition  II  above 
10  db.  These  results  are  also  given  in  Table  II  and  we  see  that  they  are 
in  good  agreement  with  the  results  obtained  from  the  relative  number  of 
events  detected.  It  must  be  remembered  however  that  this  is  a crude  com- 
parison as  we  have  included  in  each  average  low  frequency  events  from  90° 
and  high  frequency  events  from  20°  whose  losses,  as  we  shall  see,  are 
greatly  different.  Also,  in  subsequent  sections,  we  shall  uncover  substan- 
tial biases  in  the  data.  Finally  we  are  comparing  the  infinite  velocity 
beams  to  a LASA  beam  set  which  has,  as  of  this  writing,  been  substantially 
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improved  by  the  addition  of  an  infinite  velocity  subarray  beam,  together  with 
a redistribution  of  the  outer  5 subarray  beams  to  better  cover  short  wave- 
lengths. Thus  the  losses  in  Table  II  are  less  than  would  be  observed  in  a 
comparison  with  this  improved  beam  set. 

Much  ci  the  change  in  threshold  is,  rf  course,  due  to  the  changes  in 
noise  reduction  as  the  number  of  eleme  its  in  each  subarray  changes.  To 
evaluate  the  effects  of  signal  loss  as  a function  of  frequency  on  discrimina- 
tion using,  e.g.,  cepstral  analysis  or  spectral  ratios  (0.4-0. 8 Hz/l.4-1,8  Hz), 
it  is  necessary  to  allow  for  this  noise  loss  reduction.  In  the  following 
section  we  shall  discuss  the  signal  loss  as  a function  of  frequency  from  a 
theoretical  point  of  view;  and  show  that,  when  combined  with  the  expected 
noise  reduction,  and  with  allowance  made  for  effects  of  bias,  the  predicted 
S/N  loss  is  in  agreement  with  the  results  given  in  Table  II,  and  with  plots 
of  the  data  presented  in  other  ways.  Thus  we  will  be  able  to  give  theoretical 
estimates  of  signal  and  noise  loss  which  are  compatible  with  the  data  and 
which  enable  us  to  make  estimates  of  the  potential  losses  in  detection  and 
discrimination  which  will  result  from  the  formation  of  unphased  sums  for 
subarray  beams . 


-7- 


- 


THEORETICAL  SUBARRAY  SIGNAL  RESPONSE 


In  Figures  I and  2 we  see  that  the  array  response  out  to  K = K = 0.143 
(7.0  km  wavelength,  1/K  = A = v/f)  of  the  AO-16  and  4 element  subarrays.  We 
see  that  the  response  is  almost  perfectly  circular,  a fact  which  is  emphasized 
by  ligure  3 which  shows  the  mean  and  standard  deviatior  of  the  array  responses 
tor  4,  7,  10,  13  and  16  elements  averaged  over  azimuth.  Formulas  for  the  array 
response  are  given  in  the  Appendix.  Also  in  Figure  3 the  dashed  lines  give 
the  loss  at  1 Hz  relative  to  a 15.5  km/sec  beam  directed  to  the  proper  azimuth. 

In  Figure  4 we  see  the  phase  response  for  the  unphased  sums  (infinite 
velocity  subarray  beam)  for  each  experiment  as  a function  of  signal  wavelength 
and  azimuth  of  approach.  Only  angles  from  0-60°  are  plotted  since  the  response 
should  be  periodic  because  of  the  hexagonal  symmetry  of  all  the  subarray 
designs. 

As  a sidelight  it  is  of  interest  that  all  the  amplitudes  responses  are 
similar  n normalized  by  the  maximum  intersensor  spacing.  This  result  is 
shown  in  Figure  5. 

Figures  6,  7 and  8 give  the  mean  signal  loss  of  the  unphased  sums  as  a 
function  of  array  size  and  period  for  the  event  distances  30°,  60°,  and  90° 
respectively.  We  see  that  for  the  periods  plotted  the  signal  distortion 
would  not  be  great  in  the  AI  at  a distance  of  90°,  but  that  a distance  of  30°, 
any  subarray  larger  than  4 elements  would  strongly  distort  the  spectrum  and 
raise  the  threshold  of  any  discriminant  requiring  high-frequency  information, 
e.g.  a cepstral  depth  estimate,  by  0.2-0. 3 more  magnitide  units. 
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experimental  results  as  a function  of  wavelength,  period  and  velocity 

Figures  9 through  14  give  the  signal  to  noise  loss  as  a function  of  wave- 
length for  experiments  1-6  respectively.  The  wavelength  was  calculated  by  the 
formula  A = v/f  = vT  where  v is  determined  by  standard  tables  from  the  distance 
to  the  event,  and  where  T is  the  reported  period  for  the  event.  Superimposed 
on  these  figures  is  the  theoretical  mean  array  response  from  Figure  3,  dis- 
placed downwards  so  that  il  passes  through  the  median  point  in  the  interval 
12-18  km.  This  wavelength  interval  was  chosen  because  the  Partition  I sub- 
array beams  were  set  at  a velocity  of  15.5  km/sec  (A  = 15.5  km  at  1.0  Hz), 
and  thus  the  loss  of  the  unphased  sums  at  this  velocity  should  be  the  true 
loss  relative  to  a perfectly  formed  subarray  beam.  This  will  be  in  error  by 
less  than  1 db  because  of  the  slight  average  error  in  azimuth  for  the  closest 
subarray  beam.  (This  median  loss  is  tabulated  in  column  4 of  Table  III.) 

The  vertical  shift  of  the  solid  curves  in  Figures  9-14  then  represents  the 
estimated  relative  noise  reduction  for  this  array,  and  the  'urves  then  gives 

the  S/N  loss  (instead  of  the  signal  loss)  expected  at  all  wavelengths  relative 
to  a perfectly  formed  beam. 

Before  a comparison  can  be  made  between  theory  and  observations,  however 
allowance  must  be  made  for  the  fact  that  the  subarray  beams  In  Partition  I ar’e 
not  perfectly  in  phase  for  very  low  and  very  high  velocity  events.  We  have, 
therefore,  shown  a second  line  in  FiBures  9-14  which  corrects  for  the  db 
difference  of  the  loss  due  to  the  unphased  sums  (infinite  velocity  beam)  and 
that  of  the  15.5  km/sec  beam.  This  difference  curve  is  tangent  to  the 
unphased  sum  (infinite  velocity)  loss  curve  at  15.5  km/sec  since  the  15.5 
km/sec  beam  has  zero  loss  there,  and  the  difference  at  other  wavelengths  is 
also  equal  to  the  difference  in  Figure  3. 

We  see  that  for  10-16  elements  per  subarray  (Experiments  4-6,  Figures  12-14) 
the  differential  response  line  falls  below  the  data  points  for  7 < A < 10. 

This  is  as  would  be  expected  due  to  the  substantial  loss  in  S/N  for  low- 
wavelength  events  on  infinite  velocity  beams  which  results  in  their  not  being 
detected  at  all  in  Partition  II,  and  therefore  results  in  their  not  being 
Plotted  at  all  on  the  Figures.  Thus  we  see  that  the  data  as  presented  have 
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a substantial  bias  at  low  wavelengths,  and  simple  extrapolation  or  interpre- 
tations could  lead  to  substantial  errors  in  judgment. 

For  4-7  sensors  per  subarray  (Experiments  2,  3;  Figures  10,  11)  the 
ditterential  response  matches  the  data  points  satisfactorily;  suggesting  that 
the  small  differential  loss  of  about  4 db  leads  to  little  bias. 

This  line  of  argument  leads  one  to  expect  that  the  data  points  would 
cluster  about  the  differential  line  for  1 sensor  per  subarray  (Experiment  1, 
figure  9)  instead  of  falling  below  it  as  observed.  The  only  explanation  we 
can  give  for  this  behavior  is  that  the  reflection  from  the  free  surface 
arriving  at  the  500  foot  deep  seismometer  is  more  out  of  phase  with  the  direct 
arrival  than  it  is  at  the  200  foot  deep  seismometers.  The  effect  will  be 
largest  for  the  high  frequency  events  at  large  incidence  angles  (low 
velocities);  i.e.,  t-hose  from  close  in. 

Figure  9 also  offers  an  opportunity  to  evaluate  the  effects  of  signal 
contamination  of  the  noise,  the  possibility  of  which  is  discussed  in  the 
INTRODUCTION.  The  effect  should  be  greatest  in  experiment  1 where  the 
difference  in  S/N  ratio  between  partitions  is  also  the  greatest.  Or;  Figure  9 
we  have  indicated  that  the  median  value  of  S/N  loss  for  all  events,  for  those 
for  which  the  S/N  in  Partition  I was  less  than  20  db,  greater  than  20  db,  and 
greater  than  32  db  was  5.6,  5.2,  5.8  and  5.8  db  respectively.  Thus  it  would 
appear  that  with  respect  to  the  overall  average,  bias  from  noise  contamination 
is  _<  0.2  db,  a neg.ligable  effect. 

As  we  have  seen,  however,  the  bias  due  to  non— detection  at  low  wave- 
lengths has  a major  effect  in  the  data  for  10,  13  and  16  elements  per  subarray, 
and  this  fact  is  reflected  in  Figures  15  and  16  for  10  and  16  elements  per 
subarray  respectively  where  as  a result  of  the  bias  the  ful]  range  of  the 
variation  of  S/N  loss  with  period  cannot  be  seen.  Still,  it  is  clear  in  the 
case  of  16  elements  per  subarray  that  for  10°  < A < 50°  there  is  a substantial 
increase  in  signal  loss  with  frequency,  and  that  this  increase  is  in  substan- 
tial agreement  with  the  losses  pi'p  icted  from  transformations  of  the  dashed 
line  in  Figure  14,  when  allowance  is  made  for  the  effects  of  bias. 
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Similarly,  in  Figure  17  we  see  that  regression  lines  are  in  good  agree- 
ment with  the  theoretical  lines  from  Figure  9-14  evaluated  at  1.0  Hz  except 
at  short  distances  where  the  effects  of  bias  become  important. 

Figures  18-20  for  A = 30°,  60°,  and  90°  respectively  give  the  expected 
S/N  loss  relative  to  a perfectly  formed  suDarray  beam  for  0.5,  1.0,  2.0  Hz 
as  a function  of  number  of  sensors  per  subarray.  These  Figures,  derived  from 
the  solid  lines  in  Figures  9-14  correspond  to  the  signal  loss  Figures  6-8. 

An  example  of  the  use  of  these  Figures  would  be  to  consider  the  use  of  a 
short-period  discriminant  at  a distance  of  60°.  From  Figure  8 we  see  that 
for  a 10-element  array  a short  period  0.6  second  spectral  amplitude  would 
have  to  be  corrected  upward  by  4 db  before  computation  of  a spectral  ratio; 
while  by  reference  to  Figure  20  we  see  that  the  threshold  magnitude  for 
application  of  the  discriminant  would  have  to  be  increased  by  about  7.5  db 
or  .37  magnitude  units.  It  is  worth  noting  that  the  S/N  loss  relative  to 
the  16-element  array  for  0.5  Hz  is  certainly  overestimated  because  the  noise 
correlation  is  higher  at  0.5  H<_  than  in  the  band  0.8-2. 5 Hz,  thus  giving  the 
larger  subarrays  less  of  an  advantage  over  the  small. 

The  losses  derived  from  Figures  6-8  and  18-20  are  the  ones  to  consider 
for  planning  purposes,  since  the  new  subarray  beamset  in  LBS  151  now  imple- 
mented at  SDAC  has  6 subarray  beams  including  one  unphased  sum  (infinite 
velocity).  The  full  teleseismic  velocity  space  out  to  7 km/sec  is  covered 
with  a maximum  loss  of  3 db  at  1 Hz. 

In  Table  HI  we  give  a careful  comparison  of  our  observed  results  with 
theory.  In  column  3 we  give  the  theoretical  signal  loss  for  X = 15.5  from 
Figure  3.  In  column  4 we  give  the  observed  S/N  loss  for  A = 15.5  from  the 
db  value  at  X = 15.5  of  the  lines  in  Figures  9-14.  The  difference  in  these 
two  columns  is  the  apparent  noise  loss  in  the  frequency  band  0.8-2. 5 Hz  and 
is  given  in  column  5. 

The  theoretical  noise  loss  in  the  band  0.8— 2.0  Hz  may  be  found  by  the 
techniques  discussed  by  Blandford  and  Clark  (1971).  Briefly,  they  used  the 
observed  correlation  as  a function  of  distance  to  compute  i5,  the  average 
inter-sensor  noise  correlation  for  the  array.  Then  the  db  noise  reduction 


ib  Klven  by  10  logio  (N/(l+[N-l]p).  This  technique  applied  to  the  subarray 
Jo  igr.s  in  this  report  gives  the  theoretical  noise  reduction  lli  column  6. 
Since  the  noise  spectrum  at  LASA  is  heavily  dominated  by  low  equ.encies  this 
result  should  be  a good  estimate  of  the  noi;  e reduction  in  the  passed 

°,8“2*5  i!2*  ThuS*  subtracting  column  6 from  column  5 yields  the  theoretical 
.ninus  apparent  noise  loss  column  7.  We  see  that  the  smaller  the  array  the 
better  it  performs  in  comparison  with  theory;  although  the  differences  are 
only  significant  with  respect  to  practice  for  the  1-sensor  subarray. 

°f  the  3,9  db  discrePancy  for  experiment  1 we  can,  as  discussed  above, 
explain  only  0.2  db  by  appeal  to  (S/N)  bias.  Since  the  center  element  is 
buried  500  instead  of  200  feet  we  may  expect  that  it's  noise  level  is  less 
than  for  the  average  other  sensor.  In  fact,  examination  of  Table  22  from 
Chiburis  and  Hartenberger  (1966)  shows  that  the  average  center  sensor  is 
quieter  by  1.3  db  than  the  average  other  sensor.  This  leaves  2.4  db  to  be 
accounted  for.  We  feel  that  the  probable  explanation  for  the  remaining 
discrepancy  is  that  the  travel-time  residuals,  computed  either  by  examination 
of  the  single  center-sensor  trace  or  by  cross  correlation  of  the  subarray 
beams  (which  are  heavily  dominated  by  the  closely  spaced  sensors  around  the 
center  element  and  which  would  be  expected,  no  matter  what  the  central  con- 
centration, to  give  an  average  residual  close  to  the  residual  for  the 
geometric  center  of  the  array)  would  be  most  representative  of  the  array 
center  sensors.  The  advantage  would  be  expected  to  decrease  as  the  size  of 
the  subarray  increases;  and  this  is,  in  fact,  observed  in  Table  HI. 
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SUMMARY  AND  RECOMMENDATIONS 

Figure  6-8,  and  18-20  seem  to  be  confirmed  by  comparison  of  their 
coraiiaries  with  observation.  They  may  therefore  be  used  to  deduce  the  loss 
of  signal,  and  signal-to-noise  ratio  as  a function  of  frequency  which  would 
result  from  the  formation  of  unphased  sums  in  place  of  the  present  LASA 
beams  which  contains  6 subarray  beams  for  each  subarray. 

It  might  be  suggested  that  it  would  be  desirable  to  lose  less  than  0.1  nr 
Unlt  d6teCtion  threshold  in  the  AI.  Then  by  Figure  21  for  A = 90°  T = 1 0 
one  would  want  to  use  10  elements  in  each  beam.  The  threshold  for ’short-period 
spectral  ratio  (T  - 0.5)  discriminant  would  then  be  raised  by  about  0.2 
magnitude  units  in  the  AI.  At  A = 60°  (Kamchatka,  Novaya  Zemlya,  Tahiti) 
the  corresponding  numbers  would  be  0.15  and  0.4  magnitude  units;  while  at 
30°  (Cuba,  Guatemala,  Bering  Straits,  Amchitka)  we  would  have  0.2  and  0.5 
magnitude  units.  These  seem  to  be  relatively  reasonable  results;  there  is 
no  obvious  break  in  any  of  the  curves  which  would  lead  to  an  optimum  via 
simple  reasoning.  Seven  sensors  per  array  also  seems  reasonable;  but  since 
it  is  easier  to  cut  a sensor  out  of  the  beam  than  to  add  it,  one  might  well 
be  biased  toward  the  higher  number. 
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THEORETICAL  AMPLITUDE  LOSS.  dB 


th? T /Jeoretical  array  response  for  infinite  velocity  beam  of 

the  inner  4-16  elements  of  LASA  subarray  AO  averaged  over  azimuth 
he  mean  and  standard  deviation  are  indicated.  Also  shown  by  dashed 
lines  is  the  mean  db  difference  in  response  between  a 15.5  km/sec 
earn  directed  at  the  event  azimuth  and  an  infinite  velocity  beam 
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SENSORS  PER  SUBARRAY 
5 KILOMETER  HIAMETER 


igure  5.  Array  response  for  infinite  velocity  beams  for  the  inner 
16  elements  of  LASA  subarray  AO  averaged  over  azimuth  as  a functi 
'f  normalized  wavelength,  A/D. 


NUMBER  OF  SENSORS  PER  SUBARRAY 


Figure  6. 
of  sensors 


Theoretical  signal  loss  for  unphased  sums  as  a function  of  numb 
per  subarray  for  several  periods  and  for  A = 30°. 
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Theoretical  signal  loss  for  unphased  sums  as  a function  of  number 
per  subarray  for  several  periods  and  for  A = 90°. 
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Figure  9.  Signal  to-noise  ratio  loss  in  db  as  a function  of  apparent 
signal  wavelength  for  experiment  1.  The  solid  line  is  the  (constant) 
irray  response  of  a single  element  shifted  down  until  it  passes  through 
the  median  db  point  in  the  wavelength  interval  12-18  km.  The  dashed 
line  shows  the  response  relative  to  a 15.5  km/sec  beam  directed  at  a 
1 Hz  signal. 


Figure  10.  Signal-to-noise  ratio  loss  in  db  as  a function  of  apparent 
signal  wavelength  for  experiment  2.  The  solid  line  is  the  array 
response  of  the  inner  4 elements  shifted  down  until  it  passes  through 
the  median  db  point  in  the  wavelength  interval  12-18  km.  The  dashed 

line  shows  the  response  relative  to  a 15.5  km/sec  beam  directed  at  a 
1 Hz  signal. 
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Figure  11.  Signal-to-nois  - ratio  loss  in  db  as  a function  of  apparent 
signal  wavelength  for  experiment  3.  The  solid  line  is  the  array 
response  of  the  inner  7 elements  shifted  down  until  it  passes  through 
the  median  db  point  m the  wavelength  interval  12-18  km.  The  dashed 

line  shows  the  response  relative  to  a 15.5  km/sec  beam  directed  at  a 
1 Hz  signal. 
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Figure  12.  Signal-to-noise  ratio  loss  in  db  as  a function  of  apparent 
signal  wavelength  for  experiment  4.  The  solid  line  is  the  array 
response  of  the  inner  10  elements  shifted  down  until  it  passes  through 
the  median  db  point  in  the  wavelength  interval  12-18  km.  The  dashed 
line  shows  the  response  relative  to  a 15.5  km/sec  beam  directed  at  a 
1 Hz  signal. 


APPARENT  SIGNAL  WAVELENGTH.  KILOMETERS.  X-  v/f 


Figure  13.  Signal-to-noise  ratio  loss  in  db  as  i i unction  of  apparent 
signal  wavelength  for  experiment  5.  I'he  solid  line  is  the  array 
response  for  the  inner  13  elements  shifted  down  until  it  passes  through 
the  median  db  point  in  the  wavelength  interval  12-18  km.  The  dashed 
line  shows  the  response  relative  to  a ’.5.5  km/sec  beam  directed  at  a 
1 Hz  signal. 
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igure  14.  Signal-to-noise  ratio  loss  in  db  as  a function  of  apparent 
signai  wavelength  for  experiment  6.  The  solid  line  is  the  array 

thePmnSH  °£2he  Jnner  16  elements  shifted  down  until  it  passes  through 
1 j 6^ian  * point  in  the  wavelength  interval  12-18  km.  The  dashed8 
lTz  sJgnal  reSp°nse  relative  to  a 15.5  km/sec  beam  directed  at  a 
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figure  15.  Signal-to-noise  loss  in  db  as  a function  of  signal  period  for 
distance  intervals  for  experiment  4;  10  sensors  per  subarray.  Line  has 
been  transformed  from  dashed  line  in  Figure  12  by  assuming  loss  at  period 
T is  equal  to  loss  at  A = vT  where  v is  the  wave  velocity  appropriate  to 
the  distance. 
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Figure  17.  Plots  of  the  signal-to-noise  ratio  losses  for  LASA  beams 
containing  unphased  subarray  sums  as  compared  to  those  made  up  of 
phased  subarray  sums  for  each  event  versus  epicentral  distance  for 
each  experiment.  Both  regression  lines  and  transformed  dashed  lines 
from  Figures  9-14  with  f = 1.0  have  been  drawn  through  the  data. 
Transformation  effected  by  assuming  loss  at  distance  A is  equal  to 
loss  at  A = vT,  T = 1.0,  v as  appropriate  according  to  distance. 
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Fip.ure  18.  Signal- to-noise  ratio  loss  for  A = 30°  as  a function  of 
number  of  sensors  per  subarray  for  several  frequencies.  Taken  from 
transformations  of  the  dashed  lines  in  Figures  9-14. 
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Figure  19.  Signal-to-noise  ratio  loss  for  A = 60°  as  a function  of 
number  of  sensors  per  subarray  for  several  frequencies.  Taken  from 
transformations  of  tiie  dashed  lines  in  Figures  9-14. 
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Figure  20.  Signal-to-noise  ratio  loss  for  A = 90°  as  a function  of 
number  of  sensors  per  subarray  for  severa]  frequencies.  Taken  from 
transformations  of  the  dashed  lines  in  Figures  9-14, 
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APPENDIX 


DEVELOPMENT  OF  ARRAY  LOSS  EQUATION 


Consider  a sinusoidal  seismic  signal  of  frequency  f,  AeJ"t,  arriving 
at  the  center  of  the  subarray  at  time  t,  thus  the  response  of  the  i'th 
seismometer  is 

ju)  (t  +At  ) 

. o i 


where: 


v 

v = apparent  velocity,  and 

d = distance  vector  from  the  center  seismometer  along  the 
direction  of  wave  propagation. 

Therefore  the  response  of  the  seismometer  relative  to  the  center 

seismometer  at  t + At.  is  : 
o i 


R = 


jw(t  +At  ) 
Ae  ox 


jujAt 

Ae 


i 


and  the  average  response  of  the  unphased  subarray  sum  relative  to  the 
phased  subarray  sum  is: 


R 


N jojAt 

l e 

i= 


where, 


N = the  number  of  sensors. 


1 v 

Then  R = — l (cos  uiAt^  + jsin  wAt_^) 

x N N 

= — ( £ cos  ooAt . + j J sin  coAt.). 
i=l  1 i=l  X 


-AI-1- 


Let 


Amplitude  - { (~  J^cos  wAt.)2  + (i  J sin  wAt.)2)1/2 
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A - j - wavelength,  then; 
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and  phase  = tan 
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lo  express  the  amplitude  response  in  terms  of  db  loss  we  have 
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